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INVESTIGATION  AT  HIGH  SUBSONIC  SPEEDS  OF  A 
BODY-CONTOURING  METHOD  FOR  ALLEVIATING 
THE  ADVERSE  INTERFERENCE  AT  THE  ROOT 
OF  A  SWEPTBACK  WING1 

By  John  B.  McDevitt  and  William  M.  Haire 

SUMMARY 


A  body-contouring  method  for  alleviating  at  subsonic  speeds  the 
adverse  Interference  at  the  root  of  a  high -aspect -ratio  sweptback  wing 
was  investigated.  Several  bodies  used  in  combination  with  a  35°  swept¬ 
back  wing  were  modified,  as  proposed  by  Kiichemann,  by  shaping  the  body 
so  as  to  counteract  the  distorting  velocities  at  the  wing -body  junction. 

Although  modifying  the  body  shape  did  not  significantly  affect  the 
aerodynamic  characteristics  at  subcritical  speeds,  beneficial  results 
were  obtained  at  free-stream  Mach  numbers  above  the  critical.  Improved 
aerodynamic  characteristics  were  evidenced  by  large  reductions  of  drag, 
an  increase  in  lift-curve  slope,  and  a  reduced  change  of  pitching -moment- 
curve  slope  with  increasing  Mach  number. 

Additional  tests,  which  involved  changes  in  wing  position  relative 
to  the  body,  indicated  that  wing  position  had  little  effect  on  the  aero¬ 
dynamic  characteristics  of  the  unmodified  wing-body  combinations. 


INTRODUCTION 


The  effect  of  sweep  for  delaying  the  onset  of  compressibility  drag 
has  generally  been  somewhat  less  beneficial  than  indicated  by  simple- 
sweep  theory.  This  is  caused,  at  least  in  part,  by  an  adverse  velocity 
distribution  at  the  root  of  the  swept  wing.  The  subsonic  wing  theories 
of  references  1  and  2  indicate  that  the  velocity  distribution  at  the  root 
of  a  swept  wing  is,  because  of  reflection  at  the  plane  of  symmetry,  dis¬ 
torted  in  such  a  manner  that  the  chordwise  position  of  maximum  velocity 


1 Supersedes  NACA  RM  A5^A22  by  John  B.  McDevitt  and  William  M.  Haire, 
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is  displaced  rearward.  Furthermore,  the  maximum  velocity  at  the  root  of 
the  swept  wing  may  he  considerably  higher  than  that  for  an  infinite  yawed 
wing  and,  consequently,  the  premature  formation  of  strong  local  shocks 
in  the  vicinity  of  the  root  can  be  expected. 

The  velocity  distortion  can  be  alleviated  by  altering  either  the 
wing  geometry  near  the  root  or  by  contouring  the  body  in  the  vicinity  of 
the  wing-body  junction.  Shaping  the  sides  of  the  body  to  conform  to  the 
general  shape  of  a  streamline  on  a  yawed  wing  has  been  suggested  in 
reference  3*  Experimental  investigations  described  in  references  4  and  5 
were  based  on  this  method  and  some  beneficial  results  were  obtained.  A 
modified  swept-wing-root  profile  was  also  investigated  in  the  research 
reported  in  reference  4. 

The  design  of  wing-fuselage  junctions  for  subsonic  speeds  has  been 
given  an  exhaustive  treatment  by  Kuchemann  (refs.  6  and  7)»  In  one  of 
the  methods  proposed  by  Kuchemann,  the  body  is  represented  by  a  cylinder 
on  which  ring  vortices  are  distributed  so  that  the  induced  axial  compo¬ 
nent  of  velocity  cancels  the  swept-wing  interference  velocity.  By  inte¬ 
gration  of  the  induced  lateral  velocity,  the  radial  modification  neces¬ 
sary  to  shape  the  wing-body  junction  is  determined.  For  further  details 
concerning  the  use  of  ring  vortices  see  reference  8. 

The  present  report  gives  the  results  of  an  experimental  investiga¬ 
tion  of  this  method  at  high  subsonic  speeds.  A  summary  of  the  theoreti¬ 
cal  method  for  determining  the  required  body  modifications  and  a  detailed 
description  of  the  application  of  this  method  as  used  in  the  present 
investigation  are  presented  in  the  Appendix  of  this  report.  Although  the 
present  investigation  was  conducted  primarily  to  evaluate  a  particular 
body  contouring  method,  several  related  studies  were  conducted.  In  order 
to  obtain  information  concerning  the  effect  of  the  swept-wing  position 
relative  to  the  body  on  the  aerodynamic  characteristics,  a  body  was 
tested  with  the  wing  in  forward  and  in  rearward  positions.  A  body  was 
also  tested  with  the  wing  at  the  center  line  of  the  body  and  in  a  raised 
position.  Additional  information  concerning  the  effect  of  mutual  inter¬ 
ference  between  body-  and  wing-induced  velocity  fields  was  obtained  by 
testing  the  wing  in  combination  with  a  body  having  a  cylindrical  mid¬ 
section. 

It  should  be  recognized  that  the  present  method  of  body  contouring 
is  primarily  concerned  with  eliminating  the  interference  at  the  root  of 
a  swept  wing.  The  present  method  should  not  be  confused  with  the  methods 
of  references  8  and  9,  which  are  based  on  altering  the  axial  distribu¬ 
tion  of  cross -sectioned  area  to  minimize  the  wave  drag  at  Mach  numbers 
of  1  or  greater. 
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NOTATION 


A 

b 

CD 


-m 


aspect  ratio 
wing  span 


drag 

lift 

qS 


drag  coefficient, 
lift  coefficient, 
pitching-moment  coefficient, 
local  wing  chord 

wing  mean  aerodynamic  chord. 


pitching  moment  about  c/4 
qSc 


fb/2  2, 

4  c  ^ 

fb/2  c  dy 


chord  at  wing-body  junction  (chord  through  the  point  of  intersec¬ 
tion  of  basic  body  and  midchord  line  of  the  swept  wing) 


dQ  maximum  body  diameter 

Z  body  length  (distance  from  nose  to  theoretical  point  of  closure) 
M  free -stream  Mach  number 

Mcr  wing  critical  Mach  number  based  on  simple-sweep  concepts 

Mdes  design  Mach  number 

P  pressure  coefficient, 

(local  static  pressure)  -  (free-stream  static  pre s sure ) 

- 

1  2 

q  free-stream  dynamic  pressure,  -fVQ 


r  body  radius 

r0  maximum  body  radius 

r j  radius  of  basic  body  at  intersection  with  midchord  line  of  swept 

wing 


k 
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S  ving  area 

—  ving  thickness-to-chord  ratio 

u  streamvise  perturbation  velocity 

V0  free- stream  velocity 

x  distance  behind  body  nose 

y  lateral  distance  from  model  center  line 

Z  ordinate  of  ving  surface,  dimensionless  vith  respect  to  the  local 

ving  chord 

z  vertical  distance  above  model  center  line 

a  angle  of  attack 

p  «/T  -  M2 

%es  ^  ^  ^des2 

A  ving  taper  ratio 

A  angle  of  sweep,  positive  vhen  svept  back 

Ae  effective  angle  of  sweep,  tan  A^  =  -  ft 

pdes 

5  distance  behind  the  leading  edge  of  the  wing-body  .junction,  dimen¬ 
sionless  with  respect  to  the  ving  chord  at  the  ving-body  junction 

P  free -stream  mass  density 

Ar  radial  modification  associated  vith  the  vortex  cylinder 

Ary  modification  to  basic  body  radius  in  horizontal  plane  of  symmetry 

Arg  modification  to  basic  body  radius  in  vertical  plane  of  symmetry 
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APPARATUS  AMD  MODELS 
Apparatus 


The  tests  were  conducted  in  the  Ames  l6-foot  high-speed  wind  tunnel 
with  the  models  mounted  on  a  sting  support  as  shown  in  figure  1.  The 
normal  and  chord  forces  and  the  pitching  moment  were  measured  with  elec¬ 
trical  strain  gages  enclosed  within  the  model.  Multiple-tube  mercury 
manometers,  connected  to  pressure  orifices  in  the  model  by  flexible 
tubing,  were  photographed  to  provide  records  of  the  pressure  distribution 
on  the  model  surface. 


Models 


In  this  report,  for  reasons  of  brevity,  those  models  designed  with¬ 
out  consideration  of  the  wing-body  interference  problem  are  designated 
as  "basic"  and  those  designed  with  consideration  of  the  interference 
problem  are  designated  as  "modified."  The  wing  used  in  combination  with 
the  various  bodies  had  an  aspect  ratio  of  6,  a  taper  ratio  of  0.5,  and 
NACA  6^ A01 5  sections  normal  to  the  50-percent-chord  line,  which  was 
swept  back  35°.  (See  fig.  2.)  For  all  the  wing-body  combinations  tested, 
the  center  lines  of  the  bodies  were  located  in  the  chord  plane  of  the 
wing,  except  for  the  model  shown  in  figure  2(c)  which  was  tested  with  the 
wing  chord  plane  at  the  body  center  line  and  also  with  it  16.1  percent  of 
the  maximum  body  diameter  above  the  body  center  line. 

Configurations  having  basic  bodies.-  The  bodies  of  revolution  shown 
in  figures  2(a),  2(b),  2(c),  and  2(d)  were  shaped  in  accordance  with  the 
Sears -Haack  formula 


r 

ro 


except  for  the  afterportion  of  the  model  shown  in  figure  2(b)  where 
straight-line  elements  were  used  as  shown. 


The  model  shown  in  figure  2(e)  had  a  body  with  a  cylindrical  mid¬ 
section  approximately  twice  the  length  of  the  wing  root  chord.  The  body 
radii  forward  and  aft  of  the  cylindrical  midsection  were  computed  using 
the  Sears -Haack  formula  for  which  the  corresponding  values  for  l/2  were 
taken  as  the  body  length  forward  and  aft  of  the  cylindrical  midsection, 
respectively,  with  r0  equal  to  the  cylinder  radius. 


All  the  bodies  were  truncated,  as  indicated  in  figure  2,  in  order 
to  permit  mounting  on  the  sting.  The  fineness  ratios  of  the  bodies  were 
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computed  using  the  theoretical  length  of  the  body  to  closure  and  the 
maximum  body  diameter. 

Configurations  having  modified  bodies . -  The  bodies  chosen  for  modi- 
fication  were  those  shown  in  figures  2(a),  2(c),  and  2(d).  These  bodies 
were  contoured  in  the  vicinity  of  the  wing-body  junction  so  as  to  mini¬ 
mize  the  interference  velocities  at  the  design  Mach  number,  0.87,  which 
was  equal  to  the  predicted  critical  Mach  number  of  the  swept  wing  (based 
on  simple-sweep  concepts). 

A  summary  of  the  method  for  calculating  the  body  modifications  is 
given  in  the  Appendix.  The  resultant  bodies  were  not  bodies  of  revolu¬ 
tion  in  the  vicinity  of  the  wing-body  junction  but  were  elliptical  in 
cross  section.  Details  of  the  contouring  are  presented  in  figure  3« 

(See  also  the  photographs  in  fig.  1.) 

The  basic  bodies  chosen  for  modification  were  also  tested  without 
the  wing  in  order  that  the  predicted  body  pressures,  used  in  the  contour¬ 
ing  calculations,  could  be  compared  with  measured  body  pressures. 


TESTS  AND  PROCEDURE 


The  models  were  investigated  through  a  Mach  number  range  from  0.50 
to  0.9^,  with  a  corresponding  Reynolds  number  variation  (based  on  the 
wing  mean  aerodynamic  chord)  from  2.6x10®  to  3.Uxl06  (fig.  4).  The  test 
data  have  been  corrected  for  the  effects  of  blockage  and  tunnel-wall 
constraint  by  the  methods  of  references  11  and  12. 

The  base  drag,  computed  from  the  base  areas  of  the  bodies  and  the 
difference  between  measured  base  pressures  and  the  free-stream  static 
pressure,  has  been  subtracted  from  the  drag  measurements.  The  drag  data 
are  believed  to  be  slightly  in  error  due  to  an  interaction  of  the  normal 
force  and  pitching  moment  on  the  chord-force  measurements  of  the  balance. 
This  interaction  was  known  to  be  small  and  inconsistent  so  no  correction 
could  be  applied. 


'  RESULTS  AND  DISCUSSION 
Force  Studies 


The  lift,  drag,  and  pitching-moment  data  are  shown  in  figure  5  for 
the  nine  configurations.  Cross  plots  of  these  data,  which  summarize  the 
drag,  lift,  and  static  longitudinal  stability  of  the  models,  are  pre¬ 
sented  in  figures  6  to  9. 
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Drag.-  It  is  apparent  from  the  basic  data  shown  in  figure  5  that 
the  drag  for  the  configurations  having  modified  bodies  was  much  less 
than  that  for  the  basic  configurations  at  Mach  numbers  above  the  criti¬ 
cal  (the  predicted  critical  Mach  number  for  the  swept  wing,  based  on 
simple-sweep  concepts,  is  approximately  O.87).  Modifying  the  bodies 
also  resulted  in  reductions  in  drag  at  the  higher  lift  coefficients  for 
free-stream  Mach  numbers  slightly  less  than  O.87.  The  zero-lift  drag 
at  subcritical  Mach  numbers  was  not  greatly  affected  by  the  contouring. 

The  body  modification  increased  the  Mach  number  for  drag  diver¬ 
gence,  defined  as  the  Mach  number  for  which 


by  approximately  0.02.  (See  fig.  6.)  More  striking,  however,  was  the 
reduction  in  drag  above  the  divergence  Mach  number,  the  reduction  in 
zero-lift  drag  at  a  Mach  number  of  0.9^  being  of  the  order  of  40  percent. 

The  drag  for  the  wing  in  the  rearward  position  was  only  slightly 
less  than  that  for  the  wing  in  the  forward  position  (fig.  6(a)).  The 
configuration  having  a  body  with  cylindrical  midsection,  when  compared 
with  the  one  using  the  basic  body  of  the  same  size,  had  less  drag  at 
Mach  numbers  above  the  critical  but  slightly  greater  drag  at  subcritical 
speeds  (fig.  6(b)).  Raising  the  wing  did  not  affect  the  zero-lift  drag 
(fig.  6(c)).  The  same  effects  are  shown  for  a  lift  coefficient  of  0.3 
in  figure  7. 

Lift.-  The  variation  of  lift-curve  slope  (evaluated  at  zero  lift) 
with  Mach  number  is  shown  in  figure  8.  Although  minor  differences  in 
lift-curve  slope  for  the  various  basic  and  modified  models  are  evident 
at  low  Mach  numbers,  the  most  significant  effects  occurred  at  Mach  num¬ 
bers  near  and  above  the  critical.  Modifying  the  body  shape  resulted  in 
less  variation  of  lift-curve  slope  with  Mach  number  and  an  increase  in 
the  lift-curve  slope  at  supercritical  Mach  numbers. 

A  comparison  of  the  data  for  the  models  with  the  wing  forward  and 
with  the  wing  aft  (fig.  8(a))  indicates  that  the  lift-curve  slope  for 
the  model  with  the  wing  aft  rose  to  a  somewhat  lower  maximum  value  near 
the  critical  Mach  number  but  then  decreased  in  a  somewhat  similar  manner 
with  increasing  supercritical  Mach  number.  The  model  having  a  body  with 
a  cylindrical  midsection  had  a  slightly  lower  lift-curve  slope  at  low 
speeds  than  the  basic  body  of  the  same  size  (fig.  2(d)).  The  maximum 
slope  was  reached  at  a  Mach  number  slightly  greater  than  that  for  the 
basic  body,  but  the  rate  of  change  of  lift-curve  slope  with  increasing 
supercritical  Mach  number  was  greater  (fig.  8(b)).  Raising  the  wing  to 
a  position  above  the  center  line  of  the  body  did  not  affect  the  lift- 
curve  slope  (fig.  8(c)). 
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Moment.-  The  variation  of  the  static  longitudinal  stability- 


parameter 


’(if)  ’ 


with  Mach  number  is  presented  in  figure  9. 


CL“0 


Modifying  the  body  shape  resulted  in  smaller  changes  in  the  stability 
parameter  with  changes  in  Mach  number. 


The  longitudinal  position  of  the  wing  with  respect  to  the  body  had 
a  large  effect  on  the  stability  parameter  (fig.  9(a))*  The  change  in 
stability  occurring  with  increasing  Mach  number  was  most  pronounced  with 
the  wing  in  a  rearward  position.  The  effect  of  replacing  the  midsection 
of  the  body  by  a  cylinder  was  relatively  unimportant  (fig.  9(b)).  Rais¬ 
ing  the  wing  resulted  in  small  changes  in  the  stability  parameter 
(fig.  9(c)). 


Pressure  Studies 


The  effects  of  the  body  modification  on  the  pressure  distribution 
over  the  inner  half  of  the  wing  at  zero  angle  of  attack  are  illustrated 
in  figure  10  for  three  high  subsonic  Mach  numbers.  With  the  basic  bodies 
the  wing  pressures  near  the  root  were  less  negative  than  those  along  the 
outer  panel  of  the  wing,  but,  as  predicted  by  theory,  an  unfavorable 
distortion  of  the  pressure  distribution  occurred  near  the  root,  the  loca¬ 
tion  of  the  point  of  minimum  pressure  being  shifted  rearward. 

The  bodies  had  been  modified  in  an  attempt  to  obtain  the  same  pres¬ 
sure  distribution  near  the  root  chord  as  over  the  outer  portion  of  the 
wing.  Although  this  was  far  from  being  accomplished,  a  more  favorable 
distribution  of  pressure  over  the  root  was  obtained.  Specifically,  the 
point  of  most  negative  pressure  was  moved  from  a  position  behind  the 
root  midchord  to  a  position  forward  of  the  midchord.  Hie  influence  of 
the  body  modifications  extended  to  a  considerable  distance  from  the  root, 
and  the  absence  of  severe  pressure  gradients  near  the  trailing  edge  of 
the  inner  portion  of  the  wing  is  particularly  noteworthy. 


CONCLUDING  REMARKS 


These  tests  show  that  a  considerable  improvement  in  the  aerodynamic 
characteristics  of  wing-body  combinations,  employing  a  swept  wing  of 
large  thickness-to-chord  ratio  and  large  aspect  ratio,  can  be  obtained 
at  high-subsonic  Mach  numbers  by  contouring  the  body  in  the  vicinity  of 
the  wing-body  junction  according  to  the  ring-vortex  design  method  sug¬ 
gested  by  Kiichemann.  This  type  of  body  modification  did  not  significantly 
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affect  the  aerodynamic  charac ter i sties  at  Mach  numbers  below  the  critical 
of  the  swept  wing.  At  Mach  numbers  above  the  critical,  improved  aero¬ 
dynamic  characteristics  were  evidenced  by  large  reductions  of  drag,  an 
increase  in  lift-curve  slope,  and  a  reduced  change  of  pitching-moment- 
curve  slope  with  increasing  Mach  number.  Vertical  and  longitudinal 
changes  of  the  position  of  the  wing  on  the  body  had  little  effect  on  the 
aerodynamic  characteristics  of  the  unmodified  wing-body  combinations. 

Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  Jan.  22,  I95U 
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APPENDIX 

METHOD  OF  BODY  MODIFICATION 


The  design  objective  is  to  alter  the  body  shape  so  that  the  velocity 
distribution  at  zero  lift,  for  a  given  subsonic  Mach  number,  will  be  the 
same  at  the  wing-body  junction  as  that  on  an  infinite  yawed  wing  at  the 
same  Mach  number.  Computations  of  the  required  body  contours  are  based 
on  a  method  proposed  by  Kuchemann  (refs.  6  and  7)  in  which  linearized 
theory  is  used  throughout. 


Perturbation  Velocities  at  the  Wing-Body  Junction 


The  streamwise  perturbation  velocity  uj  in  the  plane  of  the  wing 
and  at  the  wing-body  junction  is  regarded  as  being  composed  of  the  fol¬ 
lowing  additive  components: 


vo 


(Al) 


where 

ut,  perturbation  velocity  induced  by  a  body  alone 

uw  perturbation  velocity  of  an  infinite  yawed  wing 

uc  distortion  velocity  at  the  root  of  a  swept  wing1 

um  perturbation  velocity  induced  by  modifications  of  the  basic  body 
shape 

In  order  that  the  velocity  at  the  wing-body  junction  be  the  same  as 
for  an  infinite  yawed  wing,  it  is  necessary  that  the  body  modification 
result  in  a  perturbation  velocity  of  sufficient  strength  to  cancel  the 
sum  of  the  body  perturbation  velocity  u^  and  the  wing-root  distortion 
velocity  uc$  that  is,  the  condition 


is  to  be  satisfied  along  the  wing-body  junction  (reference  chord). 
xIn  this  procedure  the  body  is  considered  to  be  a  reflection  plane. 
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Design  Mach  Number 


In  general,  the  compressibility  corrections  will  be  different  for 
the  two  velocity  ratios  uc/V0  and  %/Vq  to  be  canceled  by  means  of  the 
body  modification.  Because  the  variation  with  Mach  number  is  different 
for  these  two  velocity  ratios,  it  is  necessary  to  select  a  design  Mach 
number  for  use  in  computing  the  body  modification.  The  most  beneficial 
results  of  the  body  contouring  are  expected  at  supercritical  speeds; 
therefore,  it  seems  rational  to  use  the  critical  Mach  number  of  the  yawed 
wing  as  the  design  Mach  number,  since  this  is  the  highest  subsonic  Mach 
number  for  which  linearized  theory  can  be  expected  to  give  reasonable 
results.  The  NACA  6^2015  airfoil  used  has  a  theoretical  critical  Mach 
number  of  O.87  when  yawed  35°. 


Body-Induced  Velocities 


The  body-induced  velocities  may  be  computed  by  means  of  any  of  the 
methods  proposed  for  bodies  of  revolution  in  axially  symmetric  flow.  The 
slender-body  approximation  of  Lai tone  (ref.  13)  is  used  here.2  The  for¬ 
mula  is 


j(n) , 


(A3) 


where  S(x)  is  the  cross-sectional  area  of  the  body  and  SKW(x)  is  the 
nth  derivative  of  this  area  with  respect  to  x.  The  body  equation  for 
use  in  equation  (A3)  is,  for  the  bodies  of  this  report,  the  Sears-Haack 
formula 


r 

ro 


Kl-¥)T 


(A4) 


It  should  be  noted  that  equation  (A3)  includes  the  effect  of  com¬ 
pressibility.  Hie  first  two  terms  are  sufficient  to  give  a  reasonably 

2The  use  of  Lai tone’s  method,  which  depends  on  expanding  the  expression 
for  the  cross-sectional  area  of  the  body  in  a  power  series,  is  permis¬ 
sible  here  because  the  bodies  used  are  described  by  an  analytically 
continuous  function. 
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good  prediction  of  the  body-induced  velocities  near  the  center  section 
of  slender  bodies.  A  comparison  of  calculated  and  measured  body  pres¬ 
sures  is  presented  in  figure  11. 


Distortion  Velocity  at  the  Root  of  a  Swept  Wing 


The  velocity  near  the  root  of  a  swept  wing  is  known  to  be  distorted 
from  that  of  an  infinite  yawed  wing.  Kiichemann  has  determined  the  magni¬ 
tude  of  the  distorting  velocity  at  the  root  of  an  un tapered  wing  of 
infinite  aspect  ratio  in  incompressible  flow  to  be 


V0 


cos  A 


dZ(6) 

at 


n 


1  +  sin  A 
1  -  sin  A 


(A5) 


where  Z(i)  is  the  ordinate  of  the  airfoil  surface,  made  dimensionless 
with  respect  to  the  wing  chord,  and  A  is  the  angle  of  sweep  (positive 
when  swept  back) . 

For  compressible  flow,  the  three-dimensional  form  of  the  Prandtl- 
Glauert  transformation  (ref.  l4)  is  applied  to  equation  (A5) .  At  a 
given  design  Mach  number,  Mdes,  equation  (A5)  becomes 


.  VA.Z-1 

,  .  cos  Ae  — 

(^n.) _ sl  »  I*.-!.-* 

'  Vo  'des  *Pdes  1  -  sin  Ag 


(A6) 


where 

%es  =  /TT^7 


and  the  effective  angle  of  sweep  is  related  to  the  actual  sweep  angle 
by 


tan  Ag 


tan  A 
^des 


The  distortion  velocity  at  the  root  of  the  35°  swept  wing  used  in 
the  present  investigation  is  shown  in  figure  12.  The  effects  of  taper 
and  aspect  ratio  for  this  wing  are  believed  to  be  small  and  were  neg¬ 
lected.  Ifce  slope  of  the  airfoil  surface  was  determined  graphically. 
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Body  Modification 


The  body  is  replaced  by  a  vortex  cylinder,  that  is,  a  cylinder  upon 
which  ring  vortices  are  distributed  so  as  to  induce  an  axial -velocity 
ratio  um/vQ  of  sufficient  strength  to  cancel  the  velocity  ratios 
ub/v0  and  uc/vo.  The  vortex  cylinder  is  assumed  to  be  of  constant 
radius,  rj,  equal  to  the  body  radius  at  the  midpoint  of  the  reference 
chord  (at  the  wing-body  Junction)  and  of  length,  cj,  equal  to  the  length 
of  the  reference  chord.  (See  fig.  3*) 

Six  standard  distributions  of  vortex  rings,  suitable  for  the  design 
of  wing -body  intersection  lines,  and  the  corresponding  induced  velocities, 
are  given  in  reference  7*  The  distribution  functions  are: 

Even  functions 

~~  -  Jl  ~  (2t  -  if 

7°z 

IS-  =  t/l  -  (2i  -  l)2 

°6 

~2_  *  1  -  (2|  -  l)2 

7°8 

Odd  functions 

-T3 - (21-1)  7  1  -  (21  -  l)2 

°3 


=  -  (25  -  1)  tjf  1  -  (21  -  l)2 

7°7 


= -  (21  - 1)  r  1  -  (21  - 1)2 

7oa  L 


Here,  a  coordinate  system  is  used  in  which  (  is  made  dimensionless 
with  respect  to  the  reference  chord,  positive  when  directed  downstream 
and  with  origin  at  the  leading  edge  of  the  reference  chord.  The  func¬ 
tions  7v/70v  are  plotted  in  figure  13. 


NACA  TN  3672 


3.1+ 


The  distribution  of  ring-vortex  strength  along  the  cylinder  is 
determined  by  equating  the  induced  axial  velocity  ratio  to  the  value  of 

j  at  a  number  of  control  points  along  the  reference  chord. 

The  corresponding  radial  component  of  induced  velocity  is  then  used  to 
determine  the  body- shape  modification  Ar. 

The  axial  velocity  induced  by  the  vortex  cylinder  can  be  repre¬ 
sented  by 


■mi 


I 


(A7) 


and  the  radial  increment  Ar  by 


Ar  _  V 

^'4* 


(A8) 


where  ?y/70v>  u V*/^, ’  and  Arv*/(cj/2)  are  functions  tabulated  in 

tables  III,  IV,  and  VI,  respectively,  of  reference  7  for  various  values 
of  the  vortex-cylinder  length-to-diameter  ratio.  (Values  of  x  in  the 
tables  of  reference  7  correspond  to  values  of  2|-1  in  the  notation 
used  here.) 


Since  there  are  six  vortex  distribution  functions  to  be  used,  the 
velocity  ratio  um/vo  must  be  evaluated  at  six  chordwise  control  points. 
For  the  work  of  this  report,  control  points  at  t  =  0.1,  0.3,  0.5,  0.6, 
0.7,  and  0.9  were  chosen.  The  tables  in  reference  7,  which  are  for 


0, 


may  be  used  directly  if  the  velocity  ratio 


( 


ub  +  ^c\ 

Vo  ) 


is  reduced 


'by  3cles  and  if  a  vortex  cylinder  of  reduced  radius  3eiesrj  is  used* 

(See  ref.  14.)  Equation  (A7)  may  then,  with  the  help  of  tables  III  and  IV 
in  reference  7,  be  applied  at  each  control  point  and  the  resulting  system 
of  six  equations  can  be  solved  simultaneously  for  the  six  unknown  coef¬ 
ficients  Cy.  When  the  coefficients  Cy  are  known,  the  modification  Ar 
can  then  be  calculated  by  using  equation  (A8),  for  which  tabulated  values 
of  Arv/(cj/2)  are  to  be  found  in  table  VI  of  reference  7. 
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As  an  illustration,  the  radial  modifications  for  the  body  having  a. 
fineness  ratio  of  12  (wing  in  the  forward  position)  are  shown  in  fig¬ 
ure  14(a).  It  is  readily  seen  that  the  radial  modification  (calculated 
by  summing  the  odd  and  even  functions)  rapidly  approaches  zero  upstream 
and  downstream  of  the  reference  chord.  For  application  in  the  present 
method  of  body  contouring,  the  value  at  4  =  -0.2  was  subtracted  in 
order  to  start  the  modification  at  that  point,  and  an  arbitrary  fairing 
was  used  to  terminate  the  modification  at  4  =  1.2.  (See  fig.  14(b).) 


Additional  Modification 


In  the  preceding  development,  the  vortex  intensity  was  constant 
around  the  periphery  of  each  cylindrical  cross  section,  so  the  calcu¬ 
lated  changes  in  body  radius  resulted  in  annual  bumps,  invariant  with 
position  around  the  circle.  The  major  portion  of  the  “interference 
velocity"  um,  however,  is  attributable  to  the  wing  and  attenuates 
rapidly  in  the  vertical  direction.  An  approximation  for  the  attenuation 
at  the  top  of  the  body  was  obtained  by  considering  the  variation  with 
vertical  distance  of  the  maximum  perturbation  velocity  induced  by  a  two- 
dimensional,  unswept,  biconvex  airfoil.  According  to  linearized  theory, 
this  variation  is 


With  this  expression  as  a  guide,  the  radial  modification  at  the  top  of 
the  body  was  reduced  by  the  amount 


Arz 

Ar 


^fldesrj 

CJ 


tan-1 


CJ 

2SdesrJ 


(A10) 


and  the  contoured  body  shape  was  made  elliptical  in  cross  section  with 
the  full  radial  increment  at  the  side,  so  that  Ary  ■  Ar.  Details  of 
the  contoured  bodies  are  presented  in  figure  3. 
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A-164M 

(a)  General  view  of  swept  wing  mounted  on  a  modified  small  body  having  a 

fineness  ratio  of  9- 


A-164M 

(t)  General  view  of  swept  wing  mounted  on  a  modified  large  body  having  a 

fineness  ratio  of  9. 

figure  1.-  Photographs  of  the  models. 
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(c)  Close-up  of  intersection  of  swept  wing  with  modified  large  body  hav¬ 
ing  a  fineness  ratio  of  9- 


(d)  Close-up  of  intersection  of  swept  wing  with  basic  large  body  having 

a  fineness  ratio  of  9* 


Figure  1.-  Concluded. 


Figure  2. -Geometry  of  the  basic  models. 
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(b)  Small  body  having  a  fineness  ratio  of  9.  (Sea  fig.  2(d).) 
Figure  3  -Continued. 
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(c)  Large  body  having  a  fineness  ratio  of  9.  (See  fig.  2(c).) 
Figure  3 Concluded 
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Mach  number,  M 

Figure  4.-  Variation  of  Reynolds  number  with  Mach  number. 
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•  Satie  body,  wing  forward 
b  Satie  body,  wing  aft 

*  Modified  body,  wing  forward 
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(e)  Largo  bodjr  with  g  finnmtt  ratio  of  9. 
Figure  5  -  Concluded 
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Basic  body  (wing  forward ')- 

_1 - , - , - j— 

Basic  body  ( wing  oft) 


Modified  body  (wing  forward) 


(a)  Body  with  a  fineness  ratio  of  12 


Basic  body- 
Cylindrical  body - 


Modified  body 


(b)  Small  body  with  a  fineness  ratio  of  9. 


Basic  body— 

— - — 1 - r 

Basic  body  (raised  wing) 


Modified  body  -\Yy 


.5  .6  .7  .8  9  u 

Mach  number,  M 

(c)  Large  body  with  a  fineness  ratio  of  9. 

Figure  6 Variation  of  dreg  coefficient  at  zero  lift  with  Mach  number 
far  the  swept  wing  in  combination  with  various  bodies. 
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(a)  Body  with  a  fineness  ratio  of  !2 


(b)  Small  body  with  a  fineness  ratio  of  9. 


Mach  number,  M 

fej  Large  body  with  a  fineness  ratio  of  9. 

Figure  8  .-Variation  of  lift -curve  slope  with  Mach  number  for  the  swept 
wing  in  combination  with  various  bodies. 
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(a)  Body  having  a  fineness  ratio  of  12  (wing  forward). 
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(b)  Body  having  a  finonoss  ratio  of  12  (wing  aft). 
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(c)  Small  body  having  a  finonoss  ratio  of  9. 
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(d)  Largo  body  having  a  finonoss  ratio  of  9. 

Figaro  it-  Pressure  distribution  from  body-alone  toms  compared 
with  theoretical  pressure  distributions  calculated  by  equation 
(43)  for  0  Mach  number  of  0.82. 
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Figure  12.-  Distortion  velocity  at  the  root  of  a  35*  swept  back,  untapered 
wing  having  $4gA0/5  sections  normal  to  the  leading  edge  (Mach 
number  of  0.87). 
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(a)  Radial  modification  for  the  odd  and  even  vortex- 

distribution  functions . 


r 


(b)  Radial  modification  adjusted  to  have  zero  value  at 

-0.2. 

Figure  14.- Radial  modification  as  computed  by  the  ring - 
vortex  method  for  the  body  having  a  fineness  ratio 
of  12  (wing  in  forward  position). 
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